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Abstract 



The structure of both quasi-real and highly virtual photons is investigated using the 
reaction e + e~ — ► e + e~/i + /i~, proceeding via the exchange of two photons. The results 
are based on the complete OPAL dataset taken at e + e~ centre-of-mass energies close 
to the mass of the Z boson. The QED structure function F^ and the differential 
cross-section da/dx for quasi-real photons are obtained as functions of the fractional 
momentum x from the muon momentum which is carried by the struck muon in the 
quasi-real photon for values of Q 2 ranging from 1.5 to 400 GeV 2 . The differential cross- 
section da/dx for highly virtual photons is measured for 1.5 < Q 2 < 30 GeV 2 and 
1.5 < P 2 < 20 GeV 2 , where Q 2 and P 2 are the negative values of the four-momentum 
squared of the two photons such that Q 2 > P 2 . Based on azimuthal correlations 
the QED structure functions F^ and F^ for quasi-real photons are determined for an 
average Q 2 of 5.4 GeV 2 . 
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1 Introduction 



The investigation of the structure of the photon represents a fundamental test of the 
predictions of QED and QCD. The classical method of investigation is the measurement 
of photon structure functions in deep inelastic electron-photon scattering at e + e~ collid- 
ers. The photon couples to the electric charge and it reveals its structure in the fluctu- 
ations into virtual lepton and quark pairs. The pure QED process e + e~ — > e + e~/i + /i~, 
which mainly proceeds via the exchange of two photons, is an ideal environment free 
of QCD effects. In the phase space region under consideration the contribution of the 
exchange of Z bosons is negligible. For the largest part of the cross-section, both ex- 
changed photons are quasi-real and the electrons^ are scattered, undetected, at small 
angles. If one of the photons is highly virtual the corresponding electron is usually scat- 
tered into the acceptance of the detector and the reaction e7 — > e7*7 — > efi + fi~ can be 
described as deep inelastic electron scattering off a quasi-real photon, as illustrated in 
Figure [p. 
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e (Pi) 



^ 



e (P2) 




^ 



e stag(P2/ 



Figure 1 : Diagram of the reaction ee - 
denote the four- vectors of the particles 
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ee/x + /i . The symbols in brackets 



In this configuration the highly virtual photon, 7*, probes the structure of the quasi- 
real photon, 7, and the structure functions of the quasi-real photon can be measured. 
The differential cross-section fill 
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e*y— >e/i+/i 



dxdQ 2 



^Q 1 



[(! + (!- yf) F2(x, Q\ P 2 ) - y 2 F^(x, Q 2 , P 2 )] 
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at low values of y is sensitive mainly to the QED structure function F^- Here Q 2 and 
P 2 are the negative values of the four-momentum squared of the virtual photon and 
the quasi-real photon, respectively. The symbols x = Q 2 /2p ■ q and y = p ■ q/p\ ■ p 

1 Electrons and positrons are referred to as electrons 



Figure 2: Illustration of the azimuthal angle x f° r the reaction e + e — ► e + e 7*7 
e + e~fi + fj,~ in the 7*7 centre-of-mass system 



denote the usual dimensionless variables of deep inelastic scattering, and a is the fine 
structure constant. Due to the large statistics available at LEP, the dependence of F^ 
on the small virtuality, P 2 , of the quasi-real photon can be explored. The measurement 
of the distribution of the azimuthal angle, Xi between the electron scattering plane 
and the plane containing the muon pair in the 7*7 centre-of-mass system, as defined 
in Figure |2|, gives access to the structure functions Fj_ and Fg 0, as described in 
Section |2|. Real photons are only transversely polarised, but in the configuration where 
both photons are highly virtual and both electrons are detected, the cross-section for 
the process ee — > ee7*7* — > ee/i + /i~ receives sizable contributions from longitudinal 
photons. These contributions are large enough to be observed. 

Measurements of QED structure functions have been performed at several e + e~ 
colliders ||. Due to the clean final state, these measurements are limited mainly by 
statistics. At LEP, the structure function F% has been measured before by OPAL |l] 
using a smaller data statistics, and by DELPHI || and L3 ||. The distribution of 
the azimuthal angle was used to measure the structure function ratio \F^/F^ by 
OPAL [^] and Fl and F^ by L3 ||. The differential cross-section da/dx of the reaction 
ee — > ee/x + /i~ mediated by highly virtual photons has not been measured before. 

In the analysis presented here the full dataset of the OPAL detector taken at LEP 
in the years 1990 — 1995 at e + e~ centre-of-mass energies close to the mass of the Z 
boson is used. The structure function F^ and the differential cross-section da/dx for 
quasi-real photons are extracted in the largest kinematic range ever covered by a single 
experiment. In addition, the structure functions F]^ and F^ for quasi- real photons are 
determined. For the first time a measurement of da/dx for highly virtual photons is 
performed and the contributions of ttl and ttt to the cross-section are established. 
Here ttl and ttt are interference terms which correspond to specific helicity states of 
the photons, as described in detail in Section 0. 

The paper is organised as follows. The theoretical framework is outlined in Sec- 
tion 0. After a brief description of the OPAL detector in Section [3] the kinematics and 



event selection are detailed in Section [|, followed by the discussion of the results in 
Sections | and |6|. Conclusions are drawn in Section 0. 



2 Theoretical framework 



In this section, the formalism used to extract the differential cross-section and the 
structure functions is outlined. For the measurement of da/dx and F 2 , only the for- 
mulae integrated over the angular dependence of the p + p~ final state are relevant. On 
the other hand, the measurement of Fj_ and F^ involves the x dependence of the p^ pT 
final state. These two issues are discussed in turn. 

The general form of the differential cross-section for the reaction e(pi)e(p2) — ► 
e(pi)e(p 2 )X, which proceeds via the exchange of two photons 71(g), 72(^)5 is given by 



d 6 a = dV(ee->eeX) 

d 3 p' 1 d 3 p' 2 a 2 (q ■ p) 2 - q 2 p 2 



(pi -P2) 2 - mlm 2 e 



1/2 



E[E' 2 167r 4 q 2 p 

(4p+ + p++a TT + 2|p+-p+-|r TT cos 20 + 2p++p°V TL + 2pfp+ + a lIT 
+pTpfa LL -8\pi pi \r TL cos^ , (2) 

where X denotes a fermion anti-fermion state. Here p± and P2 represent the four- vectors 
of the incoming electrons, p[ and p' 2 the four-vectors of the scattered electrons and q 
and p the four- vectors of the exchanged photons 71(5) and 72 (p), where \q 2 \ > \p 2 \ is 
chosen. The scattered electrons have energies E[ and E' 2 , and is the angle between the 
two scattering planes of the electrons in the photon-photon centre-of-mass system. The 
cross-sections <jtt, 0tl, 0"lt and <tll and the interference terms ttt and ttl correspond 
to specific helicity states of the photons (T=transverse and L=longitudinal) [|TJ. Since a 
real photon can only have transverse polarisation, the terms where at least one photon 
has longitudinal polarisation have to vanish in the corresponding limit \q 2 \ — > or 
\p 2 \ — > 0, and these terms have the following functional form: o"lt oc q 2 , <Ttl oc p 2 , 
o"ll oc q 2 p 2 and ttl oc \/q 2 p 2 - The terms p? 1 and p? 2 , where j, k G (+, — , 0) denote the 
photon helicities, are elements of the photon density matrix which only depend on the 
four- vectors q, p, pi, p<i and on the mass of the electron, m e . They are listed in Ref. |l[. 

In the case of muon pair production, X = p + p~ , the cross-section is completely 
determined in QED. Equation |2| contains the full information, and it is sufficient to de- 
scribe the reaction in terms of cross-sections. However, most of the experimental results 
are expressed in terms of structure functions, since in the case of quark pair production 
the cross-section cannot be completely calculated and has to be parametrised by struc- 
ture functions. The relations between the cross- sect ions and the structure functions 
are given by || 



q 2 ^{q ■ p) 2 - q 2 p 2 ( 2 2 1 2 2 , 



4:7i z a q-p \ 2 
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1 fTr^rr ( T. CI 
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1 2 2\ 
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= F 2 7 -2xF 1 7 . 


- ^tl(x, q 2 ,p 2 



(3) 

In the limit where one virtuality, e.g. \p 2 \, is small and the other is large, \q 2 \ 3> \p 2 \, 
Eq. H reduces to 

where — q 2 = Q 2 , —p 2 = P 2 and z = E^/E^ is the ratio of the energy of the quasi- 
real photon to the energy of the beam electron radiating the quasi-real photon. In 
this formula the flux for quasi-real photons is expressed using the equivalent photon 
approximation, EPA M; 

A2 at „. ^^ i f-\ ~^2 i 0^,2 „-\ 

(5) 

For the experimental situation where the electron which radiates the quasi-real photon 
is not detected, the EPA is often used integrated over the invisible part of the P 2 
range. The integration boundary P min * s given by four-momentum conservation and 
P max is determined by # max , the maximum angle at which an electron carrying the 
energy E^ of the beam electrons could possibly escape detection. The integration of 
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-z) 2 1 

P 2 


2ml z' 
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pi 



the EPA formula leads to the Weizsacker- Williams, WW, approximation [10], which is 
a formula for the flux of collinear real photons: 



H/V f max d 2 7V 

_ / Ap2 U • X '. 



P 2 j2 

± max A 

dP 2 — 
dz J P 2 dzdP 2 

a 



2n 



1 + (1 Z) ln gnax_ 2 2 



p2 ~ ""c ~ I p2 p2 

""" ^ min max/ J 



(6) 



where P 2 in = -2—, and P 2 ^ = (1 - z)E 2 9 2 max . 



There are two potential problems with the approach of Eq. ^ which are avoided in the 
analysis presented here by using Eq. [| Firstly, if only one photon is highly virtual 
and one electron is detected, the use of the WW approximation is not adequate for 
the measurement of the P 2 dependence of the structure function Pr 7 for the quasi-real 
photon since then the dependences on P 2 are inconsistently treated because the P 2 
dependence of the EPA is integrated out, whereas the full dependence on F^x, Q 2 , P 2 ) 
on P 2 is kept. In this case the EPA should be used. Secondly, if both photons are 
highly virtual and both electrons are detected, even the EPA is not applicable, since it 
is valid only for small values of P 2 . 



If both photons are highly virtual, Eq. fj can be evaluated in the limit Q 2 3> m, 
and P 2 ^> m 2 , leading to 



dV 



d 3 pid 3 P2 



O" 



(q-pf 



q 2 p 2 



E[E 2 16vr 4 gV 

CTTT + 0"TL + 0"LT + 0"LL 



(Pi -P2) 2 -m 2 e m 2 e 
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1/2 



4/^ + P 2 + 



-ttt cos 20 — 4ttl cos ( 



(7) 



If the interference terms ttt and ttl are independent of 0, the integration over of 
the terms containing cos and cos 20 vanishes, and the cross-section is proportional to 
cttt + 0"tl + ^lt + o"ll- In this case, Eq. [| can be used to define an effective structure 
function of virtual photons, as performed in Ref. [II 
interference terms formally depend only on Q 2 , P 2 



The total cross-sections and 
on the invariant mass squared, 
W 2 , of the muon system and on the mass of the muon, m 2 However, there is a 
kinematical correlation between these variables and 0, which leads to the fact that in 
several kinematical regions, like the one used in the analysis presented, ttt and ttl are 
not independent of 0. Consequently, the terms proportional to cos and cos 20 do not 
vanish even when integrated over the full range in [|l^]. The resulting contributions 
can be very large, depending on the ratios Q 2 /P 2 , Q 2 /W 2 and P 2 /W 2 . Numerical 
results of this effect can be found in Section |5|. Due to the large interference terms, 
cancellations occur between the cross-section and interference terms and therefore no 
clear relation of a structure function to the cross-section terms can be found. In this 
situation the cleanest experimentally accessible measurement is the differential cross- 
section da/dx as predicted by Eq. [7|. 

The measurement of Fj_ and F^ requires the measurement of the \ distribution. 
Expressing the differential cross-section in terms which have the same angular depen- 
dence with respect to the azimuthal angles \ and and combinations thereof, the 
differential cross-section can be written using 13 structure functions as explained in 
Ref. |12[. By integrating over all angular dependences except the x dependence, and 



factoring out the structure function F 2 , only the structure function ratios Fj^/F 2 and 
^F^/F 2 remain for the deep inelastic electron photon scattering process and these 
ratios can be obtained from a measurement of the \ distribution. For deep inelastic 
electron photon scattering, the angle \ is defined as the angle between the deeply in- 
elastically scattered electron and the muon which, in the photon-photon centre-of-mass 
frame, is scattered at positive values of cos#*, as illustrated in Figure 0. To achieve 
sensitivity to the structure function Fj_ the definition of \ is different from that used 
in Ref. 0. The old definition lead to a vanishing term, proportional to cosx, when 
integrated of cos#*. With this new definition, the integration of Eq. [| of Ref. over 
cos#* in the range —1 to 1 leads to 2 : 



dcr(e7 — > e/i + /i ) 27ra 2 

dxdQ 2 dx/2vr ~ ~x~Q* 



(l + (l-y) 2 )x 



2 The definition of F]^ in Ref. f| differs by a factor —1/2 from that used here, because in Ref. m 
the angle x is defined differently and the integration over cos 9* is performed in the range to 1 only 



r 2 



1 + P (y)(FZ/F2) cos X + -e(y)(F^/E 



cos 2% 



Here p(y) = (2 - y)y/T=y/(l + (1 - y) 2 ) and e(y) = 2(1 - y)/(l + (1 - y) 2 ) as in 
Ref. 0. Both are close to unity for small values of y. Thus, Fj^/F^ and ^Fg/F^ are 
obtained from a fit to the x distribution. By measuring F^ in addition, Fj_ and Fg can 

~ = 0. 
and 



be calculated. Equation M is based on the structure functions for real photons, P 



The formulae for the structure functions F^, 



F^ and Fg are taken from Ref. [JO 



they keep the full dependence on the mass of the muon up to terms of order m JW . 
The mass dependent formulae are significantly different, especially for the structure 
functions Fj_ from the leading logarithmic approximation as, for example, listed in 
For example for Q 2 



Ref. [14 



5.4 GeV 2 and x = 0.8, F^ from Ref. pj is about 
10% higher than F^ using the leading logarithmic form of Ref. |14j| . The functions, 
which are explicitly used for the reweighting procedure explained in Section |6], have 
the following form: 
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4m 2 
1 ^ 



W 2 



The QED cross-section, Eq. 0, which keeps the full dependence on the virtualities 
of both photons, is implemented in the Monte Carlo programs Vermaseren [|T5"|, [TBI] , 
and GALUGA |TB[. In the implementations used here, all programs only 



BDK 



contain the multiperipheral diagram shown in Figure |T|. The contributions from the 
bremsstrahlung processes JTH] are small at low Q 2 and get more important as Q 2 in- 



creases. In this analysis the bremsstrahlung processes are treated as background using 
the predictions of the FERMISV [19| and the grc4f [2(J programs, and possible inter- 
ferences are neglected. The BDK program, in addition, contains radiative corrections. 
In the analysis presented here the Vermaseren program is used to generate a large size 
event sample which is fully simulated and treated like the data. The BDK program 
is used to determine the radiative corrections and the GALUGA program is used to 
calculate the individual contributions to the cross-section in Eq. 0. 
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In the analysis presented here, the measurement of the differential cross-section 
da/dx for events where one electron is observed (singly-tagged events) is compared to 
the QED prediction of Eq. ^. The interpretation in terms of the structure function F^ 
uses the relation between the cross-sections and the structure functions, Eq. ||. The 
measurement of the differential cross-section da/dx for events where both electrons 
are observed (doubly-tagged events) is compared to the QED predictions using Eq. [7|. 
The structure functions F"^ and F^ are obtained from a fit to the \ distribution using 
Eq.H. 



3 The OPAL detector 



The OPAL detector is described in detail elsewhere PT |. Here only the subdetectors 
which are most relevant for this analysis are briefly discussed. In the OPAL right- 
handed coordinate system the x-axis points towards the centre of the LEP ring, the 
y-axis points upwards and the z-axis points in the direction of the electron beam. The 
polar angle 9 and the azimuthal angle <fi are defined with respect to the z-axis and 
x-axis respectively. 

The OPAL detector has a uniform magnetic field of 0.435 T along the beam direction 
throughout the central tracking region, with electromagnetic and hadronic calorimetry 
and muon chambers outside the coil. The small-angle silicon tungsten calorimeter 
(SW) covers the region in 9 from 25 to 59 mrad at each end of the OPAL detector. 
The unobstructed acceptance of the forward detectors (FD) covers the 9 region from 
60 to 140 mrad at each end of the OPAL detector. Both ends of the OPAL detector 
are equipped with electromagnetic endcap calorimeters (EE) covering the polar angle 
range from 200 to 630 mrad on each side. Charged particles are detected by a silicon 
microvertex detector, a drift chamber vertex detector, a large volume jet chamber and 
a set of z-chambers. The resolution of the transverse momentum for charged particles 
is cr-pjpt = a/(0.02) 2 + (0.0015p t ) 2 for central tracks, where p t is in GeV, and degrades 
for higher values of |cos#|. The magnet return yoke is instrumented with streamer 
tubes for hadron calorimetry and is surrounded by several layers of muon chambers. 



4 Kinematics and data selection 



The observed electron which radiated the photon of higher virtuality is denoted as 
'tag' and, and in doubly-tagged events, the second observed electron is called 'stag'. 
The virtualities of the photons, Q 2 and P 2 with Q 2 > P 2 , are determined from the 
energies, E tag and E staug , and polar angles, 9 tag and s t ag , of the detected electrons using 
the relations Q 2 = 2 E h E tag (1 — cos# tag ) and P 2 = 2 E h E stag (1 — cos0 sta g), where E b 
denotes the energy of the beam electrons. By measuring the two electrons and the two 
muons, or by assuming that the second electron is scattered at cos9 8tas = ±1 in the 
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case of singly-tagged events, the kinematics is constrained. This constrain is used to 
improve on the calorimetric energy measurement of the electrons. Using conservation 
of energy and longitudinal momentum, the relation 



p M + M - cos # M+M - + (2E h - Ep+p - ) cos 9 S 

'stag — COS Ptag 



T-i ^M'M "" ""H'/i ' V.--—D *-\i'V /"""-stag / ln x 

-^tag = 7 7 U^J 

COS 6W - COS V t 



is derived. This formula allows to calculate the electron energies from their polar angles, 
and the energy E ll + ll -, the momentum Pn+p- and polar angle cos 9^+^- of the muon 
pair. The dimensionless variables of deep inelastic scattering are calculated using: 

For the singly-tagged events, P 2 is much smaller than Q 2 and is therefore neglected in 
the determination of x. 

The event selection requires in addition to one or two electrons of high energy 
detected in the electromagnetic calorimeters SW, FD or EE, the presence of exactly 
two charged particles with opposite charge which are not associated with the energy 
clusters of the scattered electrons. At least one of the charged particles has to be 
identified as a muon. The samples of singly-tagged events are denoted hereafter by the 
abbreviation of the calorimeter in which the scattered electron is detected, SW, FD 
or EE. The samples of doubly-tagged events are denoted DB. If a sample is further 
subdivided in Q 2 and P 2 , the subsamples are called, for example, SW1 or FD2. The 
samples are defined in Table [l|. For singly-tagged events, all events are vetoed which 
contain electromagnetic clusters with an energy larger than a certain fraction of E^ in 
the hemisphere opposite to the one containing the observed electron. For the SW, FD 
and EE samples the different background contributions lead to energy fraction cuts of 
20%, 20% and 5%, respectively. 

Electromagnetic clusters are accepted as electrons if they fulfill the following crite- 
ria: 

1. The energy of the cluster is larger than half the energy of the beam electrons. 

2. The polar angle 9 of the cluster is in the range 28 — 55 mrad (SW), 60 — 120 mrad 
(FD) or 210 — 540 mrad (EE) with respect to either beam direction, respec- 
tively. The numbers are chosen such that the electrons are well contained in the 
calorimeters. As the cross-section for the signal events falls off more steeply with 
9 than the cross-section for the bremsstrahlung processes, the upper limit of 9 for 
the EE sample is taken to be 540 mrad. This is a tighter limit than that required 
for good containment, but is applied to the EE sample to ensure a high signal to 
background ratio. 

A track is accepted as a charged particle if it satisfies the following criteria: 

1. It has at least 20 hits in the jet chamber. 
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2. The distance of the point of closest approach to the origin in the r0-plane is less 
than 1.0 cm in the r0-plane and less than 20 cm in the z-direction. 

3. The momentum is between 0.3 and 20 GeV, and the transverse momentum rel- 
ative to the ^-direction is greater than 0.1 GeV. 

4. The polar angle of the track is within the clean acceptance of the OPAL jet 
chamber which extends to | cos#| = 0.96. 

A particle is identified as a muon if it meets the following criteria: 

1. The momentum is greater than 1 GeV. 

2. The energy deposit in the electromagnetic calorimeter associated with the track 
is less than 1.5 GeV. 

3. It produces a muon signal either in the hadron calorimeter or in the muon cham- 
bers, as described in Ref. Q. 

The SW, FD, EE and DB samples have different kinematical distributions for the muons 
and electrons and different contributions from background processes. Therefore, the 
samples have to fulfill different trigger conditions and some specific cuts are applied in 
addition to those previously described. 

1. SW sample: 

In order to ensure a trigger efficiency close to 100%, which can reliably be es- 
timated from the data, it is required that at least one muon is observed in the 
region | cos 6^1 < 0.74. This restricted range in | cos# M | however is only required 
for the measurement of F^ and da/dx and not when measuring the x distribution 
to obtain F^/F^ and \F^/ F^. This choice is made because the accuracy of the 
measurement of Fj^/F^ and ^F^/F^ is limited by the statistical error, and the 
measurement only relies on event ratios and does not use absolute cross-sections. 
Therefore in order to retain the highest possible number of events all muons up 
to | cos# M | = 0.96 are accepted for the measurement of F^/F^ and ^F^/F^ . 

2. EE sample: 

A large background comes from Z — > fi + {i~ events with a high-energy photon 
radiated by one of the muons and where this photon is identified as an electron. 
As these photons tend to be close to the muons, this background is effectively 
rejected by requiring that the distance R = ^(Ar/) 2 + (A0) 2 between the mo- 
mentum vector of the electron candidate and the momentum vector of the closest 
muon is larger than unity. Here rj = — ln(tan(6 l /2)) is the pseudorapidity with re- 
spect to the z-axis and the azimuthal angle. In addition, at least one muon with 
momentum larger than 5 GeV and a minimum invariant mass squared, W 2 > 1 
GeV 2 , is required. 

3. DB sample: 

Due to the requirement of two observed electrons, the background is much re- 
duced. The isolation requirement is loosened to R > 0.5 with respect to the 
electron candidates in the EE calorimeter. 



12 



The ee — > ee/i + /i~ events are triggered with high efficiency by the large energy deposits 
of the scattered electrons in the electromagnetic calorimeters and by the muons detected 
in the tracking devices and muon chambers. The trigger efficiency for events fulfilling 
all kinematical cuts is determined from the data using triggers which are either related 
to the scattered electrons or to the muon tracks. The trigger efficiency is studied 
separately for the individual samples. For the EE and the FD samples completely 
independent triggers for electrons and muons are available and the trigger efficiencies 
are found to be 100% and 98%, respectively, and do not vary with x. For the SW sample 
the situation is more complicated because no trigger is available which is only based 
on the electromagnetic cluster produced by the electron. Therefore the acceptance in 
polar angle for at least one of the muons has to be reduced to get a reliable trigger 
estimate, especially for large values of x which correspond to low invariant masses of 
the muon system. The trigger efficiency for the SW sample is found to be 98% to 99% 
and does not vary with x in the ranges of x used for the analysis, as listed in Table |l|. 
The error on the trigger efficiencies are conservatively taken as the difference of the 
evaluated trigger efficiencies from 100%. 

Due to different detector configurations and status requirements, the available lu- 
minosities for the various samples are different. The luminosities used amount to 67.4 
pb _1 for the SW and DB samples, to 81.3 pb _1 for the FD sample, and to 129.1 pb _1 
for the EE sample. The number of events observed, together with the Monte Carlo 
predictions, are given in Table [|. The integrated luminosities of the signal events, 
simulated with the Vermaseren program, amount to 336 pb _1 for the SW and FD sam- 
ples, 919 pb _1 for DB samples and 964 pb _1 for the EE sample. It is verified that the 
cross-sections obtained with the Vermaseren |]15],|lIJ and GALUGA [18| programs agree 



with each other to within 1% for all samples. The background processes considered 
are 77* — > t + t~ and 77* — > hadrons, based on the Vermaseren and HERWIG 



generators, Z — ► t + t~, Z — > jjl /i~, Z — > e + e~ simulated with KORALZ p3| , and all 
processes with 7 or Z boson exchange containing four fermions in the final states as 
predicted by grc4f |20] and FERMISV ||19|| . The by far dominant source of background 



is the process 77* — * t + t~ in all samples. For the EE sample the reaction Z — > t + t~ is 
of similar importance as the reaction 77* — > t + t~ . All Monte Carlo events are passed 
through the GEANT simulation of the OPAL detector [24] and are subject to the same 
analysis as the OPAL data. 



5 Results for Jr£ and dcr/dx 



The differential cross-sections dcr/dx and the structure functions F% are unfolded from 
the observed distribution in x of the data by means of a regularised unfolding tech- 



nique [25]. The data from the different samples are analysed separately. For the 
unfolding of the structure function F^x, Q 2 , P 2 ), the {Q 2 ) and (P 2 ) values listed in 
Table [I] are used. The average values of (Q 2 ) and (P 2 ) for the different samples as 
predicted by the Monte Carlo agree well with the values observed in the data. The 
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value for (P 2 ) for the SW, FD and EE samples is taken from the Monte Carlo. The 
signal definition is based on the multipheripheral diagram only and the bremsstrahlung 
diagrams are treated as background. The contribution of the bremsstrahlung diagrams 
to the cross-section is less than 0.5% for the SW and FD sample and approximately 
1.6% for the EE sample. The resolution in Q 2 is determined from the signal Monte 
Carlo. The resolution increases with increasing Q 2 and ranges from about 2.3% for 
electrons detected in the SW calorimeter to about 3.2% for electrons detected in the 
EE calorimeter. 

For each of the samples, the agreement between the data and the Monte Carlo 
predictions, for signal and background events, is checked by comparing quantities re- 
constructed from the electrons and the muons. The variables used are the energies 
-^tag, -Estagj Ei, E 2 , and angles 6 tag , 9 sta , g , Qi, 2 , of the detected electrons and of the 
first (1) and second (2) muon, and the derived quantities y, Q 2 , P 2 and W. Good 
agreement of the distributions is found both in shape and normalisation for all sam- 
ples. The agreement between the data and the Monte Carlo predictions of the SW and 
FD samples is similar to the findings of Ref . |7[ . Some examples of control distributions 
for the EE sample and the DB samples are shown in Figures ^ and |j. The data with 
their statistical errors are compared to the signal Monte Carlo with the background 
added to it before the unfolding is performed, and to the Monte Carlo distributions 
after reweighting the signal events, where the weights are obtained from the unfolding 
procedure. All distributions show a good description of the data by the sum of signal 
and background Monte Carlo events. Figure |5] shows a comparison between the data 
and the Monte Carlo prediction for the DB samples before the unfolding using statisti- 
cal errors only. Shown are the measured azimuthal angle between the two scattering 
planes of the electrons in the photon-photon centre-of-mass system and the measured 
azimuthal angle, x, which, for doubly-tagged events, is defined with respect to the 
plane containing the electron which radiated the photon of higher virtuality. Both 
angular distributions of the data are well reproduceded by the Vermaseren program. 
The data exhibit a strong dependence on which is well described by the prediction 
of Eq. [| These distributions in principle give access to several other structure func- 
tions []I2[ but, due to the low statistics, no detailed analysis of these distributions has 



been performed. 

The differential cross-sections da/dx and the structure functions F 2 are unfolded 
from the observed x distribution of the data, shown in Figure |6|. The distributions 
in Figure |B| are not corrected for radiative effects and trigger inefficiencies. Because 
the predicted x distributions are close to the observed distributions, the results of the 
unfolding are also very close to the QED prediction. The method used to unfold the 
cross-section da/dx and the structure function F% closely follows the procedure applied 
to the measurement of the hadronic structure function -F^had anc ^ * s described in detail 



in Ref. [26]. There is however a very important difference between the /x + /i~ final 
state and the hadronic final state. The resolution in x of the /i + yU~ final state which, 
as determined form the signal Monte Carlo, amounts to about 0.03, is much better 
than the resolution in x for the hadronic final state due to the good measurement of 
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the invariant mass of the muon pair. Therefore a much finer binning in x with much 
reduced correlations between the bins, could be chosen here. 



'■i 



The differential cross-sections da/dx and the structure functions F%, normalised 
by the fine structure constant a, unfolded from the data, are shown in Figure [7] and 
Figure ^| for the independent samples of singly-tagged events and in Figure || and 
Figure |10| for the combined SW and FD samples. The measured values are corrected 
both for radiative effects and for the trigger inefficiencies. The radiative corrections 
are based on the ratio of the predicted cross-sections by the BDK and Vermaseren 
programs in bins of x. The radiative corrections vary with Q 2 and x and amount up 
to about ±9%. The structure function values are given at the centres of the x bins. 
Because, in a given bin in x, the average value of F2 is different from the value of F% 
at the centre of the x bin, the results has to be corrected for this bin size effect. The 
measured average value of F% in a given bin in x, as obtained from the unfolding, is 
corrected for the bin size effect by multiplying the measured value of i 7 ^ with the QED 
prediction of the ratio of F^ at the centre of the bin and the average F% in the bin. 
In general the corrections are small and the largest corrections occur at low, and high 
values of x. The corrections are below 1% in all but the lowest and highest bins in x, 
for all samples. In the lowest and highest bins in x the correction is always positive, 
and amounts up to 7% in the lowest x bins, and in the highest x bins it is around 5%, 
with the exception of the SW1 sample, where it is largest and amounts to 19%. The 
vertical error bars show both the statistical error and the full error, which is obtained 
from the quadratic sum of statistical and systematic errors. In all unfolded results, the 
statistical error is obtained from the quadratic sum of the statistical error of the data 
events and the signal Monte Carlo events. The measured values for the cross-sections 
da/dx are listed in Tables ||— 15] and the structure functions in Tables ||— 1[ For the 
cross-section measurement of the singly-tagged events, the corrected data correspond 
to the phase space defined by y < 0.5, P 2 < 1.4 GeV 2 and the Q 2 ranges listed in 
Table [I]. The full range of W is used, except for the EE sample where W 2 > 1 GeV 2 
is required. For the doubly-tagged events, the corrected data correspond to the phase 
space defined by y < 0.5 and the Q 2 and P 2 ranges listed in Table [I]. 

The systematic error receives contributions from several sources. The determina- 
tion of x and Q 2 is based on the measurement of the muon and electron momenta. 
The uncertainties in these measurements are taken into account by shifting the recon- 
structed quantities in the Monte Carlo samples according to resolution and repeating 
the unfolding. The variations performed are the following. 

1. The transverse momentum of the muons is shifted by 0.25% for tracks in the 
region |cos# M | < 0.90 and by 5% for very forward tracks, which are less well 
measured up to | cos# M | = 0.96. 

2. The polar angle and the azimuthal angle of the muons are shifted by 0.2 mrad 
for tracks in the region | cos 9^\ < 0.90 and by 1 mrad for very forward tracks up 
to | cos 9fj,\ = 0.96. 

3. The polar angle of the observed electron is shifted by 0.3 mrad, 0.7 mrad and 
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5 mrad for electrons observed in the SW, FD and EE calorimeters, respectively. 

The differences of the results based on the central values and the results obtained 
using the shifted values are added in quadrature. The systematic errors due to the 
uncertainties in the determination of trigger efficiencies and of the radiative corrections 
are also added in quadrature. The uncertainty is dominated in almost all bins by the 
statistical error. 

The predicted structure function F^ is strongly suppressed for P 2 = (P 2 ) compared 
to P 2 = 0. The measured structure functions Figures P and |10| are distinctly different 



from the predictions for P 2 = for all values of Q 2 . In general there is good agreement 
between the data and the predictions for all ranges in Q 2 , and the corresponding x 2 
probabilities P(x 2 ) are listed in Tables 0— §]. Some small differences between the data 
and the predictions can be seen for Q 2 > 10 GeV 2 . For the individual FD samples, 
FD2 and FD3, the data shows a slightly different shape for x < 0.5 than the predicted 
differential cross section, Figure |7|(d,e), and the structure function, Figure |9|(d,e), and 
for the EE sample the predicted differential cross section, Figure 0(f), and the structure 
function Figure ^|(f), are slightly higher than what is observed in the data for all x. 

The cross-sections for the samples DB1 and DB2, unfolded from the data, are listed 
in Table |5| and shown in Figure [II], together with the predictions from the Vermaseren 
and GALUGA programs. The data are well described by both Monte Carlos using 
the full cross-section from Eq. |2|. Using the GALUGA predictions, the influence of the 
non- vanishing terms proportional to cos <f) and cos 2(f) can be seen. If these terms are 
neglected, the predicted cross-section grossly overestimates the measured cross-section. 
This shows that both terms, ttt and ttl, are present, mainly at x > 0.1, and that the 
corresponding contributions to the cross-section are negative. The contribution from 
Ttl is especially very large in the specific kinematical region of the DB samples. There 
is also good agreement between the number of events predicted and observed for the 
DB3 sample but, because of its low statistics, no cross-section is evaluated for that 
sample. 



6 Results for FX and Fj 

For the combined SW and FD samples, the quantities Fj^/F^ and ^F^/F^ are obtained 
in bins of x from the azimuthal angle distributions shown in Figure [1^. The resolution 
in x is about 20 mrad over the whole range of x- The measured distributions do not 
exhibit directly the cos x and cos 2% dependences predicted by QED in Eq. [| This is 
mainly due to the loss of muons close to the beam axis in the region | cos^l > 0.96. In 
order to extract F^/F^ and ^F^/F^ , the azimuthal angle distributions are corrected 
for this and other detector effects using a bin-by-bin correction in x for ranges of x. The 
correction factor is a given bin is obtained from the distribution of the generated x f° r 
events passing the selection cuts, where each event is weighted by 1/(1 + F^/F^ cos %+ 



16 



|-Fg/ ' F] cos2x)- The ratios Fj_/ ' F] and F^/F? used in the weighting function are 
determined from the analytical QED structure functions Eq's. ||— |TT[ The weighted 
distribution is equivalent to what would be obtained from a flat x distribution. The 
corrected distributions are shown in Figure [13]. They are fitted to the following function: 

F(x) = A(l + Bcosx + Ccos2 X ), (14) 

where A is a normalisation factor, which is left free in the fit, B = F^/F^, and 
C = ^F^/F^, for e = p = 1. The correlation coefficients between the parameters B 
and C resulting from the fits are 0.04 over the whole x range and 0.26, 0.13, —0.11 and 
—0.29 for the different bins in x, as given in Table || The measured average values 
of Fj^l ' F^ an d \F-qI ^2 i n bi ns °f x i as obtained from the fit, are corrected for the 
bin size effect as before. The measured values of ^F^/F^ are in good agreement with 
the previous OPAL measurement of hF^/F^ 0. The corrected values of F^/F^ and 
\F^/F2 are listed in Table ^ and shown in Figure |TJ] compared to the QED predictions 
for Q 2 = 5.4 GeV 2 and e = p = 1. The QED predictions for the full range in x of 
FI/F2 = —0.027 and ^F^/F^ = 0.078 are in good agreement with the measured 
values, Fl/F^ = -0.036 ± 0.027 ± 0.004 and |i^/F 2 7 = 0.061 ± 0.013 ± 0.004. The 
measured values of hF^/F^ as a function of x are significantly different from zero and 
the measured shape agrees with the QED prediction, although it is not significantly 
different from a constant. 

The sources of systematic errors are 

1. Detector effects : 

The same variations as carried out for the measurement of F% and da/dx and 
described in Section [|, have been performed. 

2. Reweighting procedure: 

The reweighting procedure is tested by using a version of the TWOGEN [27] 



generator where the structure functions F£ and F^ can be chosen arbitrarily. 
By chosing one of the two structure functions to be a constant and measuring 
the other structure function, systematic uncertainties varying between 0.03% and 
3.3% for the measurement of F^/F^ and between 0.7% and 3.8% for ^F^/F^ in 
the different x ranges are found. 
3. Background: 

As the background is not subtracted from the data, the fit to the data is a 
superposition of the fit to the signal and background distributions contributing 
with their relative fractions. To take the effect of the background distribution 
into account in the estimation of the systematic error, the x distribution of the 
background is corrected, using the same procedure as for the data, and then fitted 
like the data. The parameters B and C of this fit are weighted by the relative 
contribution of the background to the data for each range in x. They constitute 
the estimate of the systematic error, which varies from 1.4% to 4.2% between the 
low and high x range. 

The strength of the x dependence varies with the scattering angle cos 9* of the muons in 
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the photon-photon centre-of-mass system. Reducing the acceptance of cos 6* enhances 
the x dependence but, to obtain a result for FJ_ and Fg which is valid for the full range 
of cos8* the measurement has to be extrapolated using the predictions of QED. The 
measurements of Ref. || are obtained in the range | cos 6**| < 0.7, and extrapolated 
to the full range in cos#*, whereas the measurement presented here is valid for the 
full angular range | cos^*! < 1. Taking this into account, the results for F^/Fr 7 and 
|Fg /Fr 7 obtained here and the measurements of Ref. || are consistent. 

From the measurements of F^/F^ and ^F^/F^ and F 7 , the values of F^ and Fg 
are obtained in the following way. The data from the SW and FD sample are combined 
and the structure function F 7 is unfolded using the same bins as those used for the 
measurement of Fj^/F^ and ^Fg/F 7 . As Eq's. §—[TI], used for the measurement of 



FI/F2 and \F^/ F21 are only valid for P 2 = 0, the measurement of F^ is corrected 
for the effect of non-zero P 2 in the data by multiplying the result of the unfolding for 
(P?) = 0.05 GeV 2 by the ratio of F 2 7 for P 2 = and F 7 for (P 2 ) = 0.05 GeV 2 as 
predicted by QED. Then Fj_ and Fg are calculated by multiplying the measured ratios 
by F2 and corrected for the bin size effect as before. The corrected values for F^, 



F^ and Fg are listed in Table [It], and F^ and Fg are shown in Figure [15]. The QED 



predictions for Q 2 = 5.4 GeV 2 and e = p = 1 nicely describe the data. 



7 Summary and conclusions 



The complete set of data collected by the OPAL experiment at centre-of-mass energies 
close to the mass of the Z boson are used to extract information about the QED 
structure of the photon. The structure function F 7 and the differential cross-section 
da/dx for quasi- real photons are measured in the Q 2 range from 1.5 to 400 GeV 2 , the 
largest kinematic range ever covered by a single experiment. The predictions of QED 
are found to be in good agreement with the data, and the predicted P 2 suppression of 
the structure function F, 7 is clearly observed. 

For the first time, the detailed quantitative analysis of the QED structure of the 
photon has been extended to highly virtual photons. Due to the non-vanishing inter- 
ference terms ttl and ttt in the kinematical region studied, strong cancellations in 
the differential cross-section occur between the cross-section and interference terms. 
Consequently, no clear relation of a structure function to the cross-section terms can 
be found. Therefore the differential cross-section for the reaction ee — > ee/i + /i~, pro- 
ceeding via the exchange of two highly virtual photons is measured instead, for the 
range 1.5 < Q 2 < 30 GeV 2 and 1.5 < P 2 < 20 GeV 2 . The QED predictions using the 
GALUGA program show good agreement with the data for the full cross-section and 
also the presence of the interference terms r TL and ttt in the data. 

The azimuthal correlations between electrons and muons are used to extract the 
structure functions Fj_ and Fg in bins of x. The results are consistent with those 
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already published and with the QED predictions. The measurements presented here 
supersede the earlier structure function results of OPAL O. [71 
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Sample 


SW1 


SW2 


SW 


Q 2 [GeV 2 ] 


1.5-3 


3-7 


1.5-7 


P 2 [GeV 2 ] 


— 


— 


— 


X 


0-0.97 


0-0.97 


0-0.97 


y 


0-0.5 


0-0.5 


0-0.5 


(Q 2 ) [GeV 2 ] 


2.2 


4.2 


3.0 


(P 2 ) [GeV 2 ] 


0.05 


0.05 


0.05 


data 


3259 ± 57 


2292 ± 48 


5551 ± 75 


signal 


3176 ± 25 


2273 ± 21 


5449 ± 33 


background 


47.9 ± 6.7 


60.6 ± 9.2 


108.6 ± 11.4 



Sample 


FD1 


FD2 


FD3 


FD 


Q 2 [GeV 2 ] 


6-10 


10-15 


15-30 


6-30 


P 2 [GeV 2 ] 


— 


— 


— 


— 


X 


0-0.97 


0-0.97 


0-0.97 


0-0.97 


y 


0-0.5 


0-0.5 


0-0.5 


0-0.5 


(Q 2 ) [GeV 2 ] 


8.4 


12.4 


21.0 


13.4 


(P 2 > [GeV 2 ] 


0.05 


0.05 


0.05 


0.05 


data 


1058 ± 33 


790 ± 28 


719 ± 27 


2567 ± 51 


signal 


988 ± 16 


803 ± 14 


738 ± 13 


2531 ± 25 


background 


41.4 ± 4.0 


39.5 ± 3.9 


38.2 ± 3.9 


119.1 ± 6.8 



Sample 


EE 


DB1 


DB2 


DB3 


Q 2 [GeV 2 ] 


70-400 


1.5-6 


6-30 


70-300 


P 2 [GeV 2 ] 


— 


1.5-6 


1.5-20 


1.5-20 


X 


0.1-0.9 


0-0.65 


0-0.75 


0-1 


y 


0-0.5 


0-0.5 


0-0.5 


0-0.5 


(Q 2 ) [GeV 2 ] 


130 


3.6 


14.0 


140 


(P 2 > [GeV 2 ] 


0.05 


2.3 


5.0 


16 


data 


163 ± 12.8 


111 ± 10.5 


116 ± 10.8 


8 ± 2.8 


signal 


161.9 ± 4.7 


85.0 ± 2.5 


102.1 ± 2.7 


7.9 ± 0.8 


background 


17.4 ± 2.7 


6.4 ± 3.1 


7.5 ± 1.4 


2.0 ± 0.7 



Table 1: The sample definitions in terms of P 2 , Q 2 , x and y, the average values (Q 2 ) 
and (P 2 ) as obtained from the data and Monte Carlo and the number of events selected 
from the data (data), together with the predictions of the signal Monte Carlo (signal) 
and background events (background). The SW and the FD samples are not statistically 
independent from the other samples, i.e. the SW sample contains all events from SW1 
and SW2 and the FD sample contains all events from FD1, FD2 and FD3 
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X 


SWl 


SW2 




SW 


da/dx [pb] 


da/dx [pb] 




da/dx [pb] 


0.00-0.10 


409.0 ± 23.6 ± 16.3 


164.4 ± 15.0 ± 25.0 


562.3 ± 28.3 ± 33.2 


0.10-0.20 


443.4 ± 19.7 ± 16.7 


217.4 ± 13.0 ± 


8.2 


663.0 ± 25.4 ± 22.3 


0.20-0.30 


423.2 ± 18.8 ± 16.2 


223.8 ± 12.3 ± 


8.2 


647.9 ± 24.4 ± 22.4 


0.30-0.40 


430.9 ± 18.3 ± 13.8 


220.9 ± 11.8 ± 


6.6 


650.3 ± 23.6 ± 21.4 


0.40-0.50 


377.7 ± 17.6 ± 10.9 


185.7 ± 10.1 ± 


4.9 


556.3 ± 21.6 ± 12.5 


0.50-0.60 


366.0 ± 18.3 ± 10.3 


184.8 ± 10.8 ± 


5.5 


542.4 ± 22.6 ± 15.0 


0.60-0.70 


380.1 ± 19.9 ± 12.4 


198.3 ± 11.1 ± 


4.9 


578.7 ± 24.3 ± 23.4 


0.70-0.80 


329.5 ± 20.6 ± 7.8 


181.0 ± 11.1 ± 


3.6 


500.9 ± 24.3 ± 12.3 


0.80-0.90 


302.5 ± 19.4 ± 7.4 


166.7 ± 11.1 ± 


3.6 


467.7 ± 23.9 ± 28.0 


0.90-0.97 


181.6 ± 18.7 ± 41.9 


124.3 ± 11.0 ± 13.5 


304.9 ± 22.3 ± 45.5 


(Q 2 ) 


2.2 GeV 2 


4.2 GeV 2 




3.0 GeV 2 


P(x 2 ) 


0.79 


0.26 




0.46 



Table 2: The measured differential cross-section da/dx for the SW samples. In addi- 
tion, the average Q 2 value for each sample, (Q 2 ), and the \ 2 probabilities, -P(x 2 ) 5 are 
listed. The sample definitions are given in Table |l|. The SW sample contains all events 
from SWl and SW2. The first error is statistical and the second systematic 
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X 


FDl 


FD 


da/dx [pb] 


da/dx [pb] 


0.00-0.10 
0.10-0.20 
0.20-0.30 
0.30-0.40 
0.40-0.50 
0.50-0.60 
0.60-0.70 
0.70-0.80 
0.80-0.90 
0.90-0.97 


35.4 ± 4.8 ± 2.8 

69.4 ± 5.6 ± 4.9 
66.2 ± 5.8 ± 4.0 
67.9 ± 5.5 ± 3.6 

71.2 ± 5.5 ± 3.8 

71.3 ± 5.4 ± 3.3 
65.3 ± 5.4 ± 2.5 
64.8 ± 5.3 ± 2.2 
72.1 ± 5.8 ± 3.5 

57.5 ± 5.6 ± 3.7 


61.2 ± 6.7 ± 4.2 
116.2 ± 8.1 ± 6.0 
109.6 ± 7.9 ± 5.0 
123.2 ± 8.1 ± 5.7 
111.5 ± 7.5 ± 4.4 

124.2 ± 7.5 ± 5.1 

134.5 ± 7.7 ± 5.7 

128.3 ± 7.4 ± 4.3 
124.2 ± 7.6 ± 4.8 

104.6 ± 7.0 ± 5.8 


(Q 2 ) 
P(x 2 ) 


8.4 GeV 2 
0.61 


13.4 GeV 2 
0.13 



X 


FD2 


FD3 


da/dx [pb] 


da/dx [pb] 


0.00-0.15 


20.2 ± 2.9 ± 1.6 


11.5 ± 2.7 ± 1.2 


0.15-0.30 


26.4 ± 3.0 ± 1.5 


21.9 ± 2.8 ± 1.5 


0.30-0.45 


27.6 ± 2.8 ± 1.7 


22.9 ± 2.9 ± 1.6 


0.45-0.60 


24.7 ± 2.5 ± 1.9 


26.9 ± 2.8 ± 1.4 


0.60-0.75 


35.7 ± 2.9 ± 1.9 


30.5 ± 2.7 ± 1.3 


0.75-0.90 


27.9 ± 2.8 ± 1.1 


28.9 ± 2.5 ± 1.1 


0.90-0.97 


24.5 ± 2.9 ± 1.0 


23.6 ± 3.0 ± 1.0 


(Q 2 ) 


12.4 GeV 2 


21.0 GeV 2 


P{x 2 ) 


0.07 


0.09 



Table 3: The differential cross-section da/dx for the FD samples. In addition, the 
average Q 2 value for each sample, (Q 2 ), and the x 2 probabilities, -P(x 2 ), are listed. 
The sample definitions are given in Table [|. The FD sample contains all events from 
FDl, FD2 and FD3. The first error is statistical and the second systematic 
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X 


EE 


da/dx [pb] 


0.1-0.4 
0.4-0.6 
0.6-0.8 
0.8-0.9 


2.78 ± 0.76 ± 0.28 
4.27 ± 0.58 ± 0.38 

5.79 ± 0.67 ± 0.39 
5.88 ± 0.68 ± 0.30 


(Q 2 ) 
P{x 2 ) 


130 GeV 2 

0.48 



Table 4: The measured differential cross-section da/dx for the EE sample. In addition, 
the average Q 2 value for the sample, (Q 2 ), and the \ 2 probability, -P(x 2 ) 5 are listed. 
The sample definition is given in Table [l|. The first error is statistical and the second 

systematic 



X 


DBl 


X 


DB2 


da/dx [pb] 


da/dx [pb] 


0.00-0.20 
0.20-0.40 
0.40-0.65 


9.77 ± 1.62 ± 0.78 

10.45 ± 1.26 ± 0.59 

4.34 ± 1.07 ± 0.21 


0.00-0.25 
0.25-0.50 
0.50-0.75 


5.26 ± 0.82 ± 0.99 
6.87 ± 0.78 ± 0.74 
2.75 ± 0.60 ± 0.22 


(Q 2 ) 

(P 2 ) 

P(x 2 ) 


3.6 GeV 2 

2.3 GeV 2 

0.26 


(Q 2 ) 

(P 2 ) 

P(x 2 ) 


14.0 GeV 2 

5.0 GeV 2 

0.32 



Table 5: The measured differential cross-section da/dx for the DB samples. In addition, 
the average Q 2 and P 2 values for each sample, (Q 2 ) and (P 2 ), and the \ 2 probabilities, 
P(x 2 ), are listed. The sample definitions are given in Table [I]. The first error is 
statistical and the second systematic 
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X 


SWl 


SW2 


SW 


F 1 


F 1 

r 2 


F 1 

r 2 


0.00-0.10 


0.115 ± 0.007 ± 0.005 


0.108 ± 0.010 ± 0.016 


0.113 ± 0.006 ± 0.007 


0.10-0.20 


0.219 ± 0.010 ± 0.008 


0.237 ± 0.014 ± 0.009 


0.230 ± 0.009 ± 0.008 


0.20-0.30 


0.282 ± 0.012 ± 0.011 


0.320 ± 0.018 ± 0.012 


0.300 ± 0.011 ± 0.010 


0.30-0.40 


0.347 ± 0.015 ± 0.011 


0.378 ± 0.020 ± 0.011 


0.363 ± 0.013 ± 0.012 


0.40-0.50 


0.356 ± 0.017 ± 0.010 


0.373 ± 0.020 ± 0.010 


0.364 ± 0.014 ± 0.008 


0.50-0.60 


0.400 ± 0.020 ± 0.011 


0.421 ± 0.025 ± 0.012 


0.409 ± 0.017 ± 0.011 


0.60-0.70 


0.483 ± 0.025 ± 0.016 


0.519 ± 0.029 ± 0.013 


0.507 ± 0.021 ± 0.020 


0.70-0.80 


0.491 ± 0.031 ± 0.012 


0.556 ± 0.034 ± 0.011 


0.516 ± 0.025 ± 0.013 


0.80-0.90 


0.532 ± 0.034 ± 0.013 


0.601 ± 0.040 ± 0.013 


0.574 ± 0.029 ± 0.034 


0.90-0.97 


0.308 ± 0.032 ± 0.071 


0.470 ± 0.041 ± 0.051 


0.397 ± 0.029 ± 0.059 


(Q 2 ) 


2.2 GeV 2 


4.2 GeV 2 


3.0 GeV 2 


P(x 2 ) 


0.79 


0.26 


0.46 



Table 6: The measured structure function F% for the SW samples. In addition, the 
average Q 2 value for each sample, (Q 2 ), and the x 2 probabilities, -P(x 2 )? are listed. 
The sample definitions are given in Table |I[ The SW sample contains all events from 
SWl and SW2. The first error is statistical and the second systematic 
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X 


FDl 


FD 


F 1 

r 2 


T?1 
r 2 


0.00-0.10 
0.10-0.20 
0.20-0.30 
0.30-0.40 
0.40-0.50 
0.50-0.60 
0.60-0.70 
0.70-0.80 
0.80-0.90 
0.90-0.97 


0.090 ± 0.012 ± 0.007 
0.271 ± 0.022 ± 0.019 
0.334 ± 0.029 ± 0.020 
0.409 ± 0.033 ± 0.022 
0.496 ± 0.038 ± 0.026 
0.563 ± 0.043 ± 0.026 
0.596 ± 0.049 ± 0.023 
0.687 ± 0.056 ± 0.023 
0.891 ± 0.072 ± 0.044 
0.761 ± 0.074 ± 0.049 


0.095 ± 0.010 ± 0.007 
0.264 ± 0.018 ± 0.014 
0.319 ± 0.023 ± 0.014 
0.428 ± 0.028 ± 0.020 
0.446 ± 0.030 ± 0.018 
0.558 ± 0.034 ± 0.023 
0.698 ± 0.040 ± 0.030 
0.770 ± 0.044 ± 0.026 
0.871 ± 0.053 ± 0.033 
0.795 ± 0.053 ± 0.044 


(Q 2 ) 
P(x 2 ) 


8.4 GeV 2 
0.61 


13.4 GeV 2 
0.13 



X 


FD2 


FD3 


F 1 

r 2 


F 1 

r 2 


0.00-0.15 


0.151 ± 0.022 ± 0.012 


0.117 ± 0.028 ± 0.012 


0.15-0.30 


0.297 ± 0.033 ± 0.017 


0.302 ± 0.039 ± 0.021 


0.30-0.45 


0.402 ± 0.041 ± 0.025 


0.403 ± 0.051 ± 0.029 


0.45-0.60 


0.434 ± 0.044 ± 0.034 


0.559 ± 0.058 ± 0.030 


0.60-0.75 


0.758 ± 0.062 ± 0.041 


0.782 ± 0.070 ± 0.034 


0.75-0.90 


0.723 ± 0.072 ± 0.028 


0.907 ± 0.080 ± 0.033 


0.90-0.97 


0.714 ± 0.085 ± 0.030 


0.802 ± 0.103 ± 0.033 


(Q 2 ) 


12.4 GeV 2 


21.0 GeV 2 


P(x 2 ) 


0.07 


0.09 



Table 7: The measured structure function F% for the FD samples. In addition, the 
average Q 2 value for each sample, (Q 2 ), and the \ 2 probabilities, -P(x 2 )> are listed. 
The sample definitions are given in Table |1|. The FD sample contains all events from 
FDl, FD2 and FD3. The first error is statistical and the second systematic 
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X 


EE 


F 1 

r 2 


0.1-0.4 


0.343 ± 0.094 ± 0.034 


0.4-0.6 


0.578 ± 0.079 ± 0.052 


0.6-0.8 


0.936 ± 0.109 ± 0.063 


0.8-0.9 


1.125 ± 0.130 ± 0.057 


(Q 2 ) 


130 GeV 2 


P(x 2 ) 


0.48 



Table 8: The measured structure function F^ for the EE sample. In addition, the 
average Q 2 value for the sample, (Q 2 ), and the \ 2 probability, P(x 2 ), are listed. The 
sample definition is given in Table fl]. The first error is statistical and the second 

systematic 



X 


nm 


\Fim 


x < 0.25 


0.176 ± 0.031 ± 0.010 


0.075 ± 0.025 ± 0.008 


0.25 - 0.50 


0.018 ± 0.028 ± 0.008 


0.099 ± 0.024 ± 0.010 


0.50 - 0.75 


-0.171 ± 0.029 ± 0.007 


0.081 ± 0.027 ± 0.011 


x > 0.75 


-0.228 ± 0.037 ± 0.014 


0.037 ± 0.033 ± 0.011 



Table 9: The measured values of the structure function ratios F^/F^ and ^F^/F^ for 
the combined SW and FD sample. The sample definitions are given in Table |l[ The 
first error is statistical and the second systematic 



X 


r 2 


n 


F 1 


x < 0.25 


0.249 ± 0.006 ± 0.008 


0.039 ± 0.007 ± 0.003 


0.029 ± 0.010 ± 0.003 


0.25 - 0.50 


0.523 ± 0.011 ± 0.014 


0.011 ± 0.016 ± 0.004 


0.101 ± 0.025 ± 0.011 


0.50 - 0.75 


0.738 ± 0.017 ± 0.019 


-0.122 ± 0.021 ± 0.006 


0.121 ± 0.041 ± 0.017 


x > 0.75 


0.871 ± 0.027 ± 0.021 


-0.201 ± 0.033 ± 0.013 


0.063 ± 0.056 ± 0.018 



Table 10: The measured values of the structure functions F%, Fj^ and F^ for the 
combined SW and FD sample. The sample definitions are given in Table [TJ. The first 
error is statistical and the second systematic 



28 



OPAL 



.s 

u 
u 

a 

£ 10 




i — i — i — i — I — i — i — i — i — I — i — i — i — r 



(a) 




30 40 

W [GeV] 




0.6 0.8 1 1.2 

E ta^b 



"1 — i — i — i — i — I — i — i — i — i — I — i — i — r 

-f- data (d) 

--■ Vermaseren. 
— Verm. unf. 



background J 




1 1.5 

minCe^) [rad] 



5 10 15 20 

minCEjEj) [GeV] 



Figure 3: Comparison of the data with the Monte Carlo predictions for the EE sample. 
Shown are (a) the measured invariant mass, W, (b) the energy of the detected electron 
normalised to the energy of the beam electrons, E tag /E h , (c) the polar angle of the 
muon which is scattered closer to the beam direction, min(6 l 1 , (9 2 ), and (d) the energy 
of the muon with lower energy, min(i? 1 , E 2 ). The points represent the data with the 
statistical error only. The dashed lines denote the distributions of the Vermaseren 
Monte Carlo with the background added to it before the unfolding is performed, the 
full lines are the same distributions after reweighting based on the unfolding procedure, 
and the hatched histograms represent the background contribution 



29 



OPAL 



a 

S 50 

u 

J2 40 

s 

W 30 



20 



10 







a 


_ i 


i i i i i i 


1 ' " 


-o 50 






(a)- 


V 






- 


a 






- 


-2 40 






— 


a 






■ 


V 








> 






i 


W 30 






— 


20 




.— '" 


_ 


10 




H 1- 


- 


n 


i 


1 I T I M-, T 1 I I n- 


? , ; 



0.6 



0.8 



a 

S 50 

u 
a 
ft 

3 40 

a 

m 

W 30 =-_ 



20 

10 





ill — r 



i i i i I i 



Estag^b 



(C) 



-kr-CT^Tfc^l J^!t. l I 1 



3=E 



s 


50 


'-H 




41 




a 






40 


c 




41 




^ 




w 


30 



10 15 20 25 30 
Q 2 [GeV 2 ] 



20 



"i — i — i — r 



"i — i — i — r 



(b) - 



E^U_ 



0.03 




e stag trad] 



"1 — i — i — i — i — I — i — i — i — i — I — i — i — i — r 

+ data (d) 

— Vermaseren 

— Verm. unf. 

I^H background 




15 

P 2 [GeV 2 ] 



20 



Figure 4: Comparison of the data with the Monte Carlo predictions for the DB 
samples. Figures (a,b) and (c,d) are for the DB1 and DB2 samples, respectively. 
Shown are (a) the energy of the second detected electron normalised to the energy of 
the beam electrons, E stag /E h , (b) the polar angle of the second detected electron, # stag , 
(c) the measured Q 2 , and (d) the measured P 2 . The points represent the data with 
the statistical error only. The dashed lines denote the distributions of the Vermaseren 
Monte Carlo with the background added to it before the unfolding is performed, the 
full lines are the same distributions after reweighting based on the unfolding procedure, 
and the hatched histograms represent the background contribution 
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Figure 5: Comparison of the data with the Monte Carlo predictions for the DB 
samples. Figures (a,b) and (c,d) are for the DB1 and DB2 samples, respectively. 
Shown are (a,c) the angle between the two scattering planes of the electrons in 
the photon-photon centre-of-mass system, and (b,d) the measured azimuthal angle, Xi 
defined with respect to the plane containing the electron which radiated the photon 
of higher virtuality. The points represent the data with the statistical error only. The 
full lines denote the distributions of the Vermaseren Monte Carlo with the background 
added to it, and the hatched histograms represent the background contribution 
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Figure 6: The measured x distributions for the various samples (a) SW1, (b) SW2, 
(c) FD1, (d) FD2, (e) FD3, (f) EE, (g) DB1 and (h) DB2. The sample definitions 
are given in Table |l[ The points represent the data with the statistical error only. 
The dashed lines denote the distributions of the Vermaseren Monte Carlo with the 
background added to it before the unfolding is performed, the full lines are the same 
distributions after reweighting based on the unfolding procedure, and the hatched 
histograms represent the background contribution 
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Figure 7: Differential cross-sections da/dx unfolded from the data for the various 
independent samples of singly-tagged events. The samples shown are (a) SW1, (b) 
SW2, (c) FD1, (d) FD2, (e) FD3 and (f) EE. The points represent the data with their 
statistical (inner error bars) and total errors (outer error bars). The full line is the 
differential cross-sections as predicted by the Vermaseren Monte Carlo. The tic marks 
at the top of the figures indicate the bin boundaries 
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Figure 8: Differential cross-sections da/dx unfolded from the data for the two com- 
bined samples of singly-tagged events. The samples shown are (a) SW and (b) FD. 
The symbols are as defined in Figure [j] 
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Figure 9: Structure functions F^ normalised by the fine structure constant a for the 
various independent samples of singly-tagged events. The ordering of sub-figures is 
the same as in Figure [7]. The points represent the data with their statistical (inner 
error bars) and total errors (outer error bars). The full line is the structure function 
F% normalised by the fine structure constant a for the (Q 2 ) and (P 2 ) values listed in 
Table [I], and the dashed line shows the structure function F^ normalised by the fine 
structure constant a for the same (Q 2 ) but for P 2 = 0. The tic marks at the top of 
the figures indicate the bin boundaries 
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Figure 10: Structure functions F^ normalised by the fine structure constant a for the 
two combined samples of singly-tagged events. The samples shown are (a) SW and (b) 
FD. The symbols are as defined in Figure 
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Figure 11: Differential cross-sections da/dx unfolded from the data, (a) for the DB1 
sample and (b) for the DB2 sample. The sample definitions are given in Table |l|. 
The points represent the data with their statistical (inner error bars) and total errors 
(outer error bars). The full line denotes the differential cross-sections as predicted by 
the Vermaseren Monte Carlo using the same bins as for the data. The additional three 
histograms represent the cross-section calculations from the GALUGA Monte Carlo 
for three different scenarios: the full cross-section (full line), the cross-section obtained 
for vanishing ttt (dot-dash) and the cross-section obtained for vanishing ttt and ttl 
(dash), as defined in Eq. ||]. The tic marks at the top of the figures indicate the bin 
boundaries 
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Figure 12: Measured azimuthal angle distributions for the combined SW and FD 
sample for different ranges in x. The points represent the data with the statistical error 
only, and the full line is the prediction of the Vermaseren Monte Carlo normalised to 
the number of events in the data 
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Figure 13: Corrected azimuthal angle distributions for the combined SW and FD 
sample for different ranges of x. The points represent the corrected data with the 
statistical error only, and the full line is a fit to the data based on Eq. II I 



37 



OPAL 



fa 



0.2 



fa 0.1 







-0.1 
-0.2 
-0.3 



(a) 



-^- data 



QED for Q 2 =5.4 GeV 2 and 8=p=l 



i i i i i i 



i . i . . i 



i .... i 



i . . i . i 



J_l I l_l I l_l L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 



fa 0.14 



T 



T 



T 




— QED for Q 2 =5.4 GeV 2 and £=p=l 




J I I I I I I I L 



J i i i i I i i i i I i i i i L 



I .... I 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x 

Figure 14: The measured structure function ratios F^/F^ and ^F^/F^ for the com- 
bined SW and FD sample. The points represent the data with their statistical (inner 
error bars) and total errors (outer error bars). The solid lines are the QED predictions 
for Q 2 = 5.4 GeV 2 and e = p = 1. x7<to f = 5.9/4 for i^/F 2 7 and x 2 /dof = 4.4/4 for 
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Figure 15: The measured structure functions F£ and F^ for the combined SW and 
FD sample. The points represent the data with their statistical (inner error bars) and 
total errors (outer error bars). The solid lines are the QED predictions for Q 2 = 5.4 
GeV 2 and e = p = 1. x 2 /dof = 3.5/4 for Fl and x7 dof = 3.8/4 for F% 
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